Nucleation of vortices by rapid thermal quench 
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We show that vortex nucleation in superfluid ^He by rapid thermal quench in the presence of 
superflow is dominated by a transverse instability of the moving normal-superfluid interface. Exact 
expressions for the instability threshold as a function of supercurrent density and the front velocity 
are found. The results are verified by numerical solution of the Ginzburg-Landau equation. 
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Formation of topological defects under a rapid quench 
is a fundamental problem of contemporary physics 
promising to shed a new light on the early stages of 
the evolution of the Universe. For homogeneous cool- 
ing a fluctuation-dominated formation mechanism has 
been suggested by Kibble and Zurek (KZ) Nor- 
mally, cooling is associated with an inhomogeneous tem- 
perature distribution accompanied by a phase separat- 
ing interface which moves through the system as tem- 
perature decreases. A generalization of the KZ scenario 
was suggested in Ref. for inhomogeneous phase tran- 
sitions in superfluids: if the front moves faster than 
the normal-superfluid interface a large supercooled re- 
gion which is left behind becomes unstable towards 
fluctuation-induced nuclei. 

Superfluid ^He offers a unique "testing ground" for 
rapid phase transitions Recent experiments where 
a rotating superfluid ■^He was locally heated well above 
the critical temperature by absorption of neutrons |^] 
revealed vortex formation under a rapid second-order 
phase transition. The TDGL analysis was applied to 
study a propagating normal-superfluid interface under 
inhomogeneous cooling ||^ and the formation of a large 
supercooled region was confirmed. The fluctuation- 
dominated mechanism may thus be responsible for cre- 
ation of initial vortex loops. It is commonly accepted 
that these initial vortex loops are further inflated by the 
superflow and give rise to a macroscopic number of large 
vortex lines filling the bulk superfluid. 

In this Letter we report a novel mechanism of vortex 
formation which overtakes growth of the initial loops ap- 
pearing in the supercooled region. We study the entire 
process of vortex formation in the presence of a super- 
flow using TDGL dynamics. We take into account the 
temperature evolution due to thermal diffusion. We find 
analytically and confirm by numerical simulations that 
the normal-superfluid interface becomes unstable with 
respect to transverse undulations in the presence of a su- 
perflow. These undulations quickly transform into large 
primary vortex loops which then separate themselves 
from the interface. Simultaneously, a very large num- 
ber of small secondary vortex/antivortex nuclei are cre- 



ated in the supercooled region by fluctuations resembling 
the conventional KZ mechanism. The primary vortex 
loops screen out the superflow in the inner region causing 
the annihilation of the secondary vortex/antivortex nu- 
clei. The number of surviving vortex loops is thus much 
smaller than what anticipated from the KZ conjecture. 

Model. - We use the TDGL model for a scalar order 
parameter tp ignoring the non-relevant complexity of the 
■^He-specific multicomponent order parameter: 



dtij = AV' + (1 - /(r, t))^ - IVIV + C(r, t). 



(1) 



Here A is the three-dimensional (3D) Laplace operator. 
Distances and time are measured in units of the coher- 
ence length ^(Too) and the characteristic time t^l (Too), 
respectively. These quantities are taken at temperature 
Too far from the heated bubble. The local temperature 
is controlled by heat diffusion and evolves as /(r, t) = 
Eo exp{—r'^/at)t~^/'^ where a is the normalized diffusion 
coefficient. £'o determines the initial temperature of the 
hot bubble T* and is proportional to the deposited en- 
ergy Eo such that - £q/ [C{T, - Too)e{Too){TT<jf'^] 
where C is the heat capacity. Since the deposited energy 
is large compared to the characteristic superfluid energy, 
we assume Eq ^ 1. The time at which the tempera- 
ture in the center of the hot bubble drops down to Tc is 



E^^'"". The Lan gcvin force C, with the correlator 
(CO = 2r/(5(r - Y')5{t - t') describes thermal fluctua- 
tions with a strength Tf at T^. 

The typical values of the Ginzburg-Landau parame- 
ters for Fermi liquids are: TQL{Taa) = tq/{1 — Too/Tc), 
C(Too) eo/(l - roo/Te)i/2, ^0 = hvF/2TTT, and 
To = Trh/STc- The diffusion constant a ^/£,o, ^ is 
the mean free path of a quasiparticle. In '^He, a is 
very large because £/^o ^ lO'^. The noise strength is 
~ Gi-i [1 - (^/T,)]-^/^ Gi = v{^)ilT, ^ 10^ is the 
Ginzburg number and i^(0) is normal density of states. 

Results. - We solved Eq. (j^) by the implicit Crank- 
Nicholson method. The integration domain was equal 
to 150^ units of Eq. (0) with 200^ mesh points. The 
boundary conditions were taken as = ikip with a 
constant k at the top and the bottom of the integration 
domain. This implies a superflow js = k\'4)\'^ along the 
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2-axis far away from the temperature bubble. The simu- 
lations were carried out on a massive parallel computer, 
the Origin 2000, at Argonne National Laboratory. 

Selected results are shown in Fig. |l]. One sees (Fig. 
|l|a-c) that without fluctuations (numerical noise only ||] ) 
the vortex rings nucleate upon the passage of the thermal 
front. Not all of the rings survive: the small ones collapse 
and only the big ones grow. Although the vortex lines 
are centered around the point of the quench, they exhibit 
a certain degree of entanglement. After a long transient 
period, most of the vortex rings reconnect and form the 
almost axisymmetric configuration. 

We find that the fluctuations have a strong effect at 
early stages: the vortices nucleate not only at the normal- 
superfluid interface, but also in the bulk of the super- 
cooled region (Fig. |lji-e). However, later on, small vor- 
tex rings in the interior collapse and only larger rings 
(primary vortices) survive and expand (Fig. Qf). 
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FIG. 1. 3D isosurface of 


= 0.4 for a = 400, Eo = 30 



and k = 0.5. (a-c) T/ — 0. Images are taken at times 
t = 36, 48, 80. (d-f), Tf = 0.002, t = 24, 48, 80. 

To elucidate the details of nucleation we considered an 
axi-symmetric version of Eq. (^ (depending on only r 
and z coordinates, A = + l/rdr + d^) for realistic ■^He 
parameters @: fc ~ 1, i?o > 1, and cr ^ 10^. The do- 
main was 500^ with lOOO'^ mesh points. We have found 
that without thermal fluctuations the vortices nucleate at 



the front of the normal-superfluid interface (black/white 
border in Fig. ^a-c) analogous to the 3D case. The initial 
instability is seen as a corrugation of the interface. The 
interface propagates towards the center, leaving the vor- 
tices behind. As thermal fluctuations are turned on, the 
vortex rings also nucleate in the bulk of the supercooled 
region (Fig. ^) resulting in the creation of the secondary 
vortex/antivortex pairs. We have found that the "pri- 
mary" vortices prevent the supercurrent from penetrat- 
ing into the region flUed with the secondary vortices. One 
sees that the primary vortices encircle the brighter spot 
in Fig. 1^ indicating a larger value of the order param- 
eter and thus a smaller value of the supercurrent. As a 
result the secondary vortices either annihilate with an- 
tivortices due to their mutual attraction or collapse due 
to the absence of the inflating superflow. Fig. || shows 
the number of vortices N"^ and antivortices vs time 
with and without fluctuations. Fluctuations initially cre- 
ate a very large number ~ 10* of vortices and antivortices 
in the bulk which then annihilate. The resulting amount 
of surviving vortices N = iV+ — is only weakly de- 
pendent on fluctuations. 




FIG. 2. Grey-coded images of \tp\ for axi-symmetric Eq. 
(0) for a = 5000, Eq = 50 and k — 0.5, black corresponds to 
\tf}\ = and white to |'i/>| = 1. Current is along the z-axis. 
Vortices are seen as black dots, (a-c) Tf — 0, images are 
shown for t = 40, 100, 200; (d-f) Tf = 0.002, for t = 30, 50, 200 
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Shown in Fig. ^ is the number of vortex rings N vs 
quench parameters and appUed current k. At small k N 
shows threshold behavior while becoming almost linear 
for larger k values. The deviations from linear a law 
appear close to the value of the critical current kc = 1 / VS 
for a homogeneous system. 
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FIG. 3. Number of vortices iV+ (filled 
symbols) and anti-vortices A*'" (open symbols) vs time for 
a — 5000, Eo = 50 and k = 0.5. Circles correspond to T/ = 
and squares to Tf = 0.002. Inset: iV = iV+ - A^" for Tf = 
(solid line) and Tf = 0.002 (dashed line) 

Stability of normal- super fluid interface.- Following 
Ref. 1^,0, we expand the local temperature 1 — / near Tc- 
Let us put X ^ rc^r where rc is the radius of the surface 
at which T = Tc or / = 1, i.e., = (3/2)at log(Wx/<)- 
A positive x is directed towards the hot region. We write 
1 — f{r,t) ~ —a{x — vt) where a — ~[df / dr]f=i — 2rcl(Tt 
is the local temperature gradient and v — (arg)^^ 
is the front velocity defined through the quench rate 
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We have v = [Zat - 2ri) /4rct. The 
front starts to move towards the center at t > t^, — t,nax/e 
and disappears at t — t^^^ when the temperature drops 
below Tc- The front velocity accelerates as the hot bub- 
ble collapses. Since the front radius rc is large compared 
to the coherence length, it can be considered flat. We 
chose the coordinate y parallel to the front. 

We transform to the frame moving with the velocity v 
and perform the scaling of variables x,y = (x,y)v, t — 
tv^ , ip = and u — v'^/a. The parameter u ^ 

(c^/imax)/ log^(tmax/i) IS of the Order 1 in the experi- 
ment 1^ at the initial time but grows rapidly as the hot 
bubble shrinks. In our numerical simulations, u ^ 1. 
Eq. (|l|) takes the form (we drop tildes in what follows) 

5tV = AV' + a.V--V'-|V'l'V'- (2) 
u 

The amplitude F of the current-carrying solution ip = 
Fexp{iky) satisfies 



0. 



(3) 



To examine the transverse stability of stationary solution 
to Eq. (^ we put ij: = {F -\- w) exp(z/cx) where the real 
and imaginary parts of the perturbation w = a + ib are 



A 
iB 



exp(A(g)t + iqy) 



(4) 



where q is the transverse undulations wavenumber and A 
is the corresponding growth rate, we obtain 



A^ + 2xB = dlA + d^A --A- iF'^A 

u 

AB + 2xA = dlB + d^B --B- F^B 

u 



where x = kg , A = X + , and x ^ x + uk"^ 



(5) 
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FIG. 4. (a) Number of vortex rings N as function of k for 
Eq = 50 and a = 5000 and (b) iV vs cr for fc = 0.4 and 
-Eo = 50. Dashed lines show fitting to prediction Eq. (p^. 

The eigenvalue A for x ~^ can be found as an expan- 
sion in X- A = xAi + x^A| -I- ... and similarly for A and 
B. In zeroth order in x one has Aq = 0,Bq — F. In the 
first order we derive Bi — and 



d^Ai + d^Ai - -Ai - BF^Ai = 2F. 



(6) 



The solution Ai = 2udxF is obtained by differentiating 
Eq. (^). In the second order to Eq. (|^) one has 

dlB2 + d.,B2 - -B2 - F^B2 = Aud^F + A^F. (7) 
u 

A zero mode of Eq. ^) is F . The adjoint function is 
B^ = F exp{x). Eq. (^ has a solution if the solvability 
condition with respect to the zero mode is fulfilled 



/oo 
dxFe''{Aud.^F + A2F) = 
-00 



(8) 



3 



After integration we obtain A2 = 2u. Returning to the 
original notations, we obtain the exact result 



A = g2(2ufc2-l) + 0(g4) 



(9) 



The instability occurs above the threshold value = 
(2u)-i or ~ a2/3/yi/3 _ cr"Mog(tmax/0 in the 
Ginzburg-Landau units. Since it is much smaller than 
the bulk critical value kc = 1/ \/3, it can be exceeded for 
a very small superflow. 

The eigenvalue A vs x a-nd u can be derived indepen- 
dently in the limit 1 assuming that A ^ x ^ l/it <C 
1. Substituting x = x — W7, where 7 determines the po- 
sition of the interface, we treat the terms containing A, x 
and x/u as perturbations for e = 1/m 0. For e = and 
7 > 1/4 Eqs. (||) possess a front solution. This solution 
should be matched with its asymptotics for a; > 0, and 
this match fixes the value of 7. As it was shown in Ref. 
0, for u ^ 00 the matching is possible for 7 ^ 1/4. 

For e = Eqs. §) have 2 zero modes: {A, B) = {F^,Q) 
and {A^B) = (0, F). The solvability conditions result in 
the characteristic equation for A: 



1 



ciA - Ac2X 



-0, 



(10) 



where the coefhcients ci^2, d are given in the forms of in- 
tegrals of F with the corresponding zero modes in the 
interval —oo<x<xq. The constant xq is determined 
from the condition d = because for x = there is 
always an exact solution to Eq. (|^) with A = 0. Substi- 
tution of the solutions for 7 ^ 1/4 yields ci ^ 2, C2 — > 1 
and the largest growthrate of transverse perturbations 



l/u- 



(11) 



Numerical solution of Eqs. (^) demonstrates an excellent 
agreement with the theoretical expression Eq. (pi]). 

Now we apply the above results to estimate the num- 
ber of nucleated vortices. The evolution of perturbations 
near the interface is given by the integral 



w 



dq exp[A((7)t + iqy] 



(12) 



In the case of thermal quench, the normal/superfluid 
front velocity m — > 00 as time increases, therefore, the 
limit of large u applies. For u —^ 00 one has A = 2\kq\—q^. 
The maximum growth rate is achieved at q — k and 
is simply k'^. Taking into account that it is the ther- 
mal noise which provides initial perturbations for the in- 
terface instability, one derives from Eqs.(ojl2) (|t«|) ^ 



^/Tf exp[A;^t + iky]. The number of vortices is estimated 
as iV = r^k, where tq is the radius of the front where the 
perturbations w become of the order of one. The time in- 
terval to corresponding to = 1 is tp ^ fc^^ log(T, 



Vortices have no time to grow if to 
then finds: 



For To one 



= (3/2)cr(tmax - to) log(tmax/(<max " ^o)) ■ (13) 

The number of vortices with logarithmic accuracy is 

N^kro^ ^El'^'^ivs/vcf - (3^ \og{TJ^)/ET (M) 

where /3 = const, while Vs and Vc are the imposed and 
critical GL superflow velocity, respectively. This estimate 
is in agreement with the results of simulations, see Fig. 
^. Eq. (^^ exhibits a slow logarithmic dependence of 
the number of vortices at the interface on the level of 
fluctuations and agrees with the results presented in Fig. 
3 §. For CT ~ 103,£'o - 10^ - 10^ which is close to 
the experimental values of the parameters. Our analysis 
results in about 10 surviving vortices per heating event. 
It is consistent with Ref. Q where as many as 6 vortices 
per neutron were detected. 

In conclusion, we have found that the rapid normal- 
superfluid transition in the presence of superflow is domi- 
nated by a transverse instability of the normal/superfluid 
interface propagating from the bulk into the normal re- 
gion. This instability produces primary vortex loops 
which then separate from the interface. Simultaneously, 
a large number of vortex/antivortcx pairs are created by 
fluctuations in the bulk of the supercooled region formed 
after the collapse of the hot bubble. The primary vortex 
loops screen out the superflow and cause annihilation of 
the vortex/antivortex pairs in the bulk. The number of 
surviving vortices is determined by superflow-dependent 
optimum wavevector of the interface instability. 
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Even for Tf = Q there is a small level of fluctuation due to 
roundoff effects in numerical integration. For single preci- 
sion it may generate noise with the strength Tj ~ 10~^. 
Eq. (^^ describes only effect of interface instability where 
as Fig. 3 displays outcome of whole nucleation process 
both in the bulk and the interface. 
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